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Abstract 

In this work we analyze the locahzation of fermions on a brane embedded in five- 
Q ■ dimensional, warped and non-warped, space-time. In both cases we use the same non- 

. linear theoretical model with a non-polynomial potential featuring a self-interacting 

scalar field whose minimum energy solution is a soliton (a kink) which can be con- 
tinuously deformed into a two- kink. Thus a single brane splits into two branes. The 
behavior of spin 1/2 fermions wavefunctions on the split brane depends on the coupling 
00 '. of fermions to the scalar field and on the geometry of the space-time. 



: 1 Introduction 
o . 

The idea that our Universe can be reahzed inside a domain wall embedded in a (4,1)- 
dimensional world [Ij has provided many creative ways to solve the hierarchy of interactions 
I problem [2]- [3] and the cosmological constant problem in worlds with large extra dimen- 
^ ' sion, without resorting to the compactification of the extra dimension. It has also been 
shown |5j that the effective gravitational potential between two particles recovers the New- 
tonian behavior, since one has localization of gravitons on a thin brane in five- dimensional 
space-time with a warped geometry and the cosmological constant is related to the brane 
tension. Localization of matter (spin- zero, spin-1/2 and spin-3/2) in the Randall- Sundrun 
framework was shown to be possible, under certain conditions over the brane tension [B]. 
The localization of spin 1/2 fermions on thin branes is due to a soliton, via the mechanism 
created by Jackiw and Rebbi [7] to demonstrate the phenomenon of fermion charge frac- 
tionization. By its turn, the domain wall which we would live in, according to the scenario 
proposed by Rubakov and Shaposhnikov \T\, is the topological defect provided by the very 
same soliton. 

The Rubakov-Shaposhnikov scenario pLj has been extended to five- dimensional warped 
space-time by means of a self-interacting scalar field, or a set of (self-)interacting scalar fields 
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also coupled to gravity [8]. In those nonlinear models, also inspired on previous studies on 
the stabilization of gravity fluctuations on domain walls in supergravity theories [9] , the thick 
branes solutions are minimum energy configurations (in many cases they are Bogomol'nyi- 
Prasad-Sommerfield (BPS) solutions) that separates the space in two patches with a peculiar 
warp-factor whose asymptotic behavior is an anti-de Sitter (AdSs) space. 

In some models the minimum energy configurations are wide solitons [10]- [12], also called 
double- wall. Double- wall configurations also appear as solutions, degenerate in energy [T3] . 
of a model with two interacting scalar fields [ll]-[15], and they have been called degenerate 
Bloch branes. The model used in [14J also comprises critical Bloch branes [13] which are 
analogues of the extreme domain walls described in [9]. Such a variety of solutions [16] 
is due to an arbitrary constant of integration of the equations of motion. The continuous 
deformation from a single brane (domain wall) to a double-brane, up to an extreme one, by 
varying the constant of integration, is similar to the phenomenon of brane splitting discussed 
in [IT], in analogy to first-order phase transitions in condensed matter systems. Such a 
transition is usually approached by using a sixth degree polynomial potential, as was done 
in [T7], which is characterized by the increase of a disordered phase (a wet phase in surface 
physics) or an increase of the thickness of the domain wall (brane). As the temperature 
of the system approaches a certain limit, one has the appearance of two interfaces between 
the disordered phase and the ordered ones, that is, the formation of a double-wall. That 
phenomenon has been called brane splitting [17]. As the temperature of the system goes 
towards a critical one, a complete domination of the disordered phase (complete wetting) 
might happen, which we interpret, in brane worlds scenario, as the formation of an extreme 
brane. Such a phase transition can also be described by effective non-polynomial potentials 
[18] constructed from the model with two scalar fields [l^, [IE], [I3] with the advantage of 
having a model whose minimum energy solutions are the BPS ones. 

The localization of fermions on a double-wall (double-brane) in warped space-time has 
been studied in [19] and in pO]. It has been shown that double- wall can localize massless 
fermions. Notwithstanding, a close inspection on the behavior of the zero-mode eingenfunc- 
tion reveals that it is peaked just in the region between the branes and has tails inside the 
branes, where the possibility of detecting massless fermions is very small. It can be noticed 
that this behavior is due to the Yukawa coupling, of fermions to the soliton. We think 

that the peak of the zero-mode eingenfunction should follows the trend of the brane, that is, 
it should also split. In this work we show that this behavior can be achieved by means of a 
convenient coupling of the fermion to the scalar field which is reminiscent from N = 1 super- 
symmetry (SUSY), namely VF00'l'\E', where W^^ is the 2nd derivative of the superpotential, 
of the field theory model, with respect to the field taken at the BPS configuration. 

We develop the calculations by using one of the models obtained in [18j, but they 
could also be done by using any other nonlinear model which admits two-kink solutions 
for the scalar field. We also analyze the localization of massless and massive fermions in the 
Rubakov-Shaposhnikov framework [Ij, that is, in a non- warped space-time. In this latter 
case one clearly sees that the coupling 0'I^\1/ provides an effective single well potential where 
the fermions would be trapped in, but that well is just in the region between the branes, 
while the coupling W^f^^p^'^ provides a double-well potential with wells in the cores of the 
branes, as it should be. In the case of fiat space-time, the Numerov method is used to obtain 
the massive localized modes, with emphasis on the SUSY inspired coupling with which there 
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is room to prepare a chiral mixed state of quasi-degenerate massless and tiny mass fermion 
states. We also show that such a mixed state can tunnel between the two branes with time of 
tunneling inversely proportional to the tiny mass which, in its turn, decreases as the distance 
between the walls increases. We would like to warn the reader that the brane splitting is not 
considered here as a dynamic process; in fact, the distance between the walls is one of the 
parameters of the nonlinear potential that could be seen as dependent on the temperature 
and the tunneling of fermions is analyzed at a fixed distance between the walls, that is, at a 
given temperature close to the critical temperature for the formation of an extreme wall. 

In the next section we show the main features of the nonlinear model coupled to gravity 
in 5-dimensional space-time [I8]. In the third section we are concerned to the localization of 
massless fermions on a split brane in warped space-time. In the fourth section we deal with 
the localization and tunneling of massive fermions in flat 5-dimensional space-time. A few 
remarks on fermion localization on double-walls are left to the last section. 



2 The model, the double-brane and the warp factor 

The model we are going to deal with includes a self-interacting scalar field coupled minimally 
to gravity with one extra space dimension, denoted by r. The action that leads to Euler- 
Lagrange for the scalar field and Einstein equations is given by 

S = J d^xdy^M (^-^-R + - V{<P)^ , (1) 

where g = Deti^gab) and the metric is 

ds^ = gahdx^dx^ = e^'^^'^^Jx^'dx" - dr\ a, 6 = 0, 4, (2) 

where t]^,^ is the Minkowski metric and e'^^'^^^ is the warp factor, which is supposed to depend 
only on the extra dimension r. The Greek indices run from to 3. 
We consider that the potential V{(j)) can be written as 



n,,.u^jmy.Uw(,)f, (3) 



2 

2\ d(j) J 3 

In this case the BPS solution of the following first-order differential equations 



are also solutions of the second-order differential equations of motion in the static limit. 
By taking W{(j)), which we call the superpotential, given by 



W{(f)) = 2/i 



3 2 



(5) 



with / = \/W—a? and b < —a < 0, we have 

dW{(f)) 



WM = 2f,i<i^' + f + b./^TP), (6) 
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and the solutions of the first-order differential equations (jl]) can be found straightforwardly. 
The classical configuration for the scalar field is given by 



cosh(2yuar) — b/ f ' 

where the (lower)upper sign stands for (anti-)kink configuration. The reference [18] has more 
details on the soliton profile, as well as on the behavior of the warp factor in this model. As 
has been noted in the second paper of reference P^, the solution ([7]) can be conveniently 
written as 

= ±— [tanh(yua(r + L)) + tanh(/ia(r — L))], (8) 



where we have identified b = — cos h(2L)/v/cosh(2L)2 - 1. The above expression can be seen 
as a merging of two solitons (two kinks), namely 0± = f [±1 + tanh(yua(r =p L))] , which are 
localized at r = —L and r = +L. 

In flat two-dimensional space-time one constructs such a scalar model as the bosonic 
sector of a = 1 SUSY model with the potential given by only the first contribution at 
the right hand side of the equation In such a circumstance one can clearly see that 
the potential has two minima at = ±a and one additional minimum appears at = as 
b = —a. This transition in the potential is reflected in the behavior of the soliton in ([8]), 
which is continuously deformed from one kink into a two-kink profile and finally in one of 
the solutions, 0±, for b = —a (in this critical case L is replaced by an arbitrary constant of 
integration). A first-order phase transition takes place. It is characterized by the growing 
of the disordered phase with the formation of two interfaces separating the disordered phase 
from the ordered ones, up to the complete wetting in the critical limit. The phase transition 
can also be seen from the behavior of the warp factor, discussed in |J8j. One can see 
that the warp factor separates the space along the extra dimension in two similar regions, 
whose asymptotic behavior is an AdSs space and, as L increases, a double-brane structure is 
triggered at a specific value of L. Such a double-brane structure is evident from the peaks 
of the Ricci scalar R = —(8 dPA/dr"^ + 20 {dA/drY) just at the points the core of the branes 
are localized (see Figure 1 below). Moreover, in the limit b = —a the branes are infinitely 
separated from each other, there being the formation of an extreme brane. Such an extreme 
case is discussed in references [13] and [21] where the brane is called a critical Block brane. 
In the next sections we focus only on the localization of massless fermions on single brane 
that splits into two brane, that we also call double-brane or double-wall. 



3 Localization of fermions: warped space-time 

In this section we analyze the localization of massless fermions on the brane described in 
the previous section by taking into account also the curvature of the space. The problem 
is approached by searching normalized solutions of the Dirac equation for fermions coupled 
with the scalar field by a general Yukawa coupling ^'F(0)\l/, where F{(j)) is a functional of 
the field 0(r) taken at the classical solution. The functional form of F{(j)), as well as the 
dependence of the metric on r, is crucial for the possible localization of massless fermions 
on the brane, as we discuss below. 
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The fermion action defined as 

S f ermion 

leads to the Dirac equation 



j d^Xy/^ {^iVDa^il - (9) 



liT^'Da - = 0. (10) 

The gamma matrices satisfy the algebra {F", F''} = 2g°''', and can be defined in the irreducible 
representation as F'^ = e~^^^'''y^, F^ = —i'j^ in terms of the 4x4 gamma matrices, 7", in 
fiat space-time. By means of the definition of the covariant derivative Da = {da + uja) and 
the metric in ([2]), one finds that the only non- vanishing components of the connection Ua 
are o;^ = |e'^(c?rA)7^75 and that the Dirac equation turns out to be written as 

{iYd^ + e^(")7^[9, + 2drA{r)] - e^(")F(0)}^(x, r) = 0. (11) 

To ease the separation of variables one can resort to the chiral decomposition (x stands for 
the four space-time coordinates) 

^(x, r) = ^ i'Ln{x)aLn{r) + ^ ipB.n{x)aRnir), (12) 

n n 

where ipL{R)n{x) are the chiral modes which satisfy •y^ipLnix) = —ipLn{x), •y^ipRnix) = ipRn{x) 
and also the four-dimensional massive Dirac equations i'y^d^ipL{R)n = '"^n''/'_R(L)n-Then, we 
find the following differential equations for the r-dependent scalar parts of the spinor 



a 



+ [2A' - F{(f))]aRn = -rrine-^aLn, (13) 



and 

a'Ln + [2A' + F(0)]ain = m„e-^ajj„, (14) 

where the prime stands for the first-derivative with respect to r. The equations above can 
be put in a more familiar form, namely 

R'^ - F{(f))Rn = -mnC-^Ln and L'^ + F{4>)Ln = mnC'^Rn, (15) 

by using the redefinitions anni'^) = e^^^*^''^i?„(r) and ai„(r) = e~^'^*^^^L„(r). 

Equations (fT5|) are equivalent to the equations for the components of a spinor describing a 
massless fermion in 1-1-1 dimensions subject to a mixing of scalar and vector potentials. The 
time- independent equation for a fermion under such potentials can be written as Hip{r) = 
Eijj{r), with the Dirac Hamiltonian given (in natural units) by if = + Vs{r)ai + K,(r), 
where p = —id/dr is the momentum operator, ai and o"2 are the two off-diagonal Pauli 
matrices, Vs{r) = —F{(f)) (note that F{(f)) is a function of r) is the scalar potential and 
K(r) = mne~^^^^ is the vector potential. Analogously, Ln{r) and Rn{r) play the role of the 
upper and lower components, for the fermion zero- mode in 1-|-1 dimensions. As a matter of 
fact, the scalar potential can be seen as a position-dependent fermion mass. As one knows, 
many examples of such systems were already solved in the literature [22], but only a few 
potentials may support bound states. 
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In the brane world scenario, by its turn, massive (m„ 7^ 0) as well as massless (m„ = 0) 
fermions might be localized inside the brane, depending on the functional coupling F{(f)). 
The issue of localization of massive modes on the brane is not going to be discussed in this 
section, even though, we would like to remark that the analysis of the possible massive spec- 
trum of bound states is not as simple as the search for localized massless modes [21], |23j . 
Moreover, normalizable ai{(2,)„(r)— functions is guaranteed if they satisfy the ortonormaliza- 
tion relations 



00 PCX) 

00 J —00 

00 

e^^^'^aLnir)aRmir)dr = 0. (16) 

-00 

From these relations one can note that the warp factor is crucial to determine the localization 
of whatever mode. 

Particularly, for the massless mode one finds that 

afioW = NR^exp[-2A{r) + J F{r')dr'] 

«LoW = NL,exp[-2A{r)- J F{r')dr'] , (17) 

where we have used F{r) = F{(f){r)) and Nr(_l) ^ are constants of normalization that are 
found according to the normalization relations 

/oo /"OO 
^-AM-f2rF(.')^.'^^ = iN^f / e-^M-2/^^(^')'^'-'cir = 1. (18) 
'OO J —00 

One can note that the normalization conditions are determined by the asymptotic behav- 
ior of the integrands in the expressions above. In most of the examples dealing with thick 
branes the factor exp[— y4(r)] goes asymptotically to infinity given the asymptotic behavior 
of the warp factor, that is exp[2A(r — )■ ±00)] 0, hence one has to choose F(0(r)) in such a 
way that the decreasing of exp[± J F{r')dr'] is faster than the increasing of exp[— 74(r)], but 
this choice does not guarantee the normalization of both chiralities, because of the different 
signs, ±, in the exponent, that is, if aRg{r) is normahzable, aLo(r) is not and vice-versa. 
Moreover, the factor exp[— 274(r)] is symmetric around the core of the brane whereas the 
kink solution is odd under the reversion of the r coordinate around the core of the brane, 
which is usually chosen at r = 0, hence one has to choose F{(f)) as an odd function on r, in 
order to have the normalizable chiral mode even in r and with a peak on the brane. 

We have analyzed the behavior of the massless chiral modes for two cases, namely: 
F(0) = ?70(r) and F{(f)) = —rjW^^, where 77 > is a coupling constant and W^^ is the second- 
derivative of the superpotential ^ with respect to taken at the two-kink configuration 
in ([7]). The first case is the simplest Yukawa coupling of fermions to a scalar real field, 
while the second one is inspired on the coupling of fermions and bosons within a = 1 



6 



SUSY model, which is also considered in the next section. This last functional coupling is 
the one which provides the correct localization on the brane. In fact, the simplest Yukawa 
coupling also provides a localized massless left-handed mode, aLQ^r), but it does not follows 
the brane splitting, that is, while the double-wall is formed, the peak of the wavefunction is 
midway between the two walls, there being a small probability density to find such a mode 
on the core of the walls themselves. On the other hand, in the case F{(f)) = —viW^^f,, one has 
(^Lo{f) with peaks on the branes, signalizing a great probability for the massless left-handed 
mode to be found just on the branes and a small probability to be found on the bulk and 
between the walls. Those behaviors can be seen from Figure 2, where the profiles of «Lo(^) 
are shown for different values of L. Figure 2 should be confronted with Figure 1. 



4 Localization of fermions: flat space-time 

In this section we adopt the usual analysis to find fermion bound states under the action of 
the scalar field whose classical configuration is given in ([7]). Particularly, we focus on the 
scenario proposed in pLj , that is, a brane (or domain wall) immersed in a five-dimensional fiat 
space-time. The action for the fermion field is given in ([9]) with gab = Vab the Mikowskian 
metric. Da = da and the irreducible form of the gamma matrices, = 7^^, is going to 
be used. The chiral decomposition (IT^ can also be used to separate the four space-time 
variables, x^, from the variable r. Now, the r-dependent functions appearing in the chiral 
decomposition obey the following equations 

"ijn - F{(j))aRn = -ninaLn, (19) 

and 

"Ln + F{(j))aLn = m^aRn. (20) 
Particularly, for the massless mode one finds 



oiR^ir) = NR,eM+ J F{r')dr'] 

/r 
F{r')dr'] . (21) 

Now, the normalization of the massless modes depends on the asymptotic behavior of 
F{r')dr' only, hence one usually has a unpaired (isolated) chiral (left-handed or right- 
handed) zero-mode, which is the main feature for having fermion number fractionization, as 
shown in |7]. 

We again have analyzed the behavior of the massless modes by setting -F(0) = ri(f)[r) and 
F(0) = —rjWtptf,, with r] > 0. In the first case we have 

a^^(r) = and aL,(r) = A^Lo(Cosh2L + Cosh2r)-''/l (22) 

As in the previous section, the function aLo('") symmetric in r and has a peak at r = 0, 
hence the massless left-handed mode is localized in the region between the branes with a 
very small probability density at the core of the branes. 
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For F{(f)) = —rjW^fjj one finds 



ajig{r) = and aLo(^) = ^Lq (sech^(r + L) + sech^(r — L 



)) 



(23) 



From Figure 3 we can note that olqI'") is symmetric in r, has no nodes and exhibits peaks at 
the cores of the branes. One can also note that the probability density to find the left-handed 
massless mode in the midway between the branes decreases as L increases. 

The adequate behavior of the massless mode is sufficient enough to consider the func- 
tional coupling F{(f)) = —rjW^^ as very convenient and has motivated us to analyze the 
consequences of such a coupling on the localization of massive modes on the split brane in 
flat space-time. As one knows, equations f lT9|) and fl20l) can be decoupled to two second-order 
differential equations, namely 



where f//j(r) = {vW^,i>y -rjW'^^ and [/^(r) = {r]W^^f + r]W^ in the case F{(f)) = -r]W^^. It 



is also known that the equations above are time-independent Schrodinger equations, whose 
corresponding Hamiltonians are superpartners of each other, that is, one has a quantum 
mechanics supersymmetry. This is formally true whatever is the functional coupling F{(j)), 
but in the case considered here such supersymmetry seems to be a reflection of a supersym- 
metry at the fundamental level. In other words, one can note that r-dependent part of the 
excitations of the scalar field {branons), picked up to quadratic terms on the fundamental 
Lagrangian density in fiat space-time (now with V^cj)) = obeys a time-dependent 

Schrodinger equation similar to the one obeyed by ai„(r) with rj = 1 , that is, with an 
effective potential given by Ueff{r) = (W^^)'^ + W^^, resulting identical mass spectrum for 
branons (bosonic excitations) and fermions . 

In Figure 4 it is shown the form of the potentials f/i?(r) and ^/^(r) for a specific value 
of L and r] = 1. For values of L close to zero, Uiir) is a single well potential, which starts 
to be deformed into a double-well potential as L approaches to a critical value Lc, that 
is determined by the condition U'l{r = 0) = 0; for L ^ Lc, dimples are observed around 
r = ±L, and a remarkable double-well is observed for the value oi L = I corresponding to 
[/^(r = 0) = 0. It worth mentioning that, although the bottom of the double-well is slightly 
raised as L increases, the width of the double-well potential also increases, signalizing the 
possible entrapment of a massive state, besides the, always present, massless one. This 
possible appearance of a massive bound state can also be seen from the behavior of Ufi{r), 
which is a single well potential whose bottom is above zero for L < I and equals to zero for 
L = I, that is, the deepness and width of Ufi{r) increase as L increases. 

One can also observe that Ui^r) = 2(2 — Ssech^r) and UR{r) = 2(2 — sech^r) for L = 
and f] = 1', such that the first potential admits two bound states and the later admits only 
one bound state. The fundamental state of t/L(r) for L = and 77 = 1 is aLo(^) — sech^r, 
whereas the first excited state is ^^^^(r) ~ sechr tanhr and the fundamental state of f/j?(r) 
for L = and rj = 1 is <yRi{r) ~ sechr. Moreover, from the expression f l23p with 77 = 1, 
one can construct an antisymmetric function as Q;Li(r) ~ sech (r + L) — sech (r — L) as an 
approximate expression for the first excited state of Uiir) when L » I. This approximation 



^Rn + UR{r)aRn = mluRn, 
-a'l^ + UL{r)aLn = rnlam , 



(24) 
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for the first excited state is commonly used to approach the discrete spectrum of double well 
potentials [21], 

We have used the results described above, together with the Numerov method to analyze 
the behavior of ai^ir), a_Ri(r) and the eigenvalue ml at intermediary values of L. Those 
behaviors are shown in Figures 5 and 6. From them one can see that a;j?^(r) is mainly 
distributed on the region between the branes, hence there is a small probability for the 
massive right-handed mode being observed inside the wells where the Universe(s) would be 
realized, whereas the probability density associated with aLi(r) is pronounced just on the 
cores of the branes. In summary, at least one massive left-handed mode is localized on 
the branes. The eigenvalue ml decreases smoothly as L increases, that is an expected result 
when one is dealing with double-well potentials in non-relativistic quantum mechanics. That 
way, one can construct a mixed left-handed massive state with both, the fundamental and 
the first excited state, which are quasi-degenerate for very large values of L, namely 

^L,n.iAt,r) = N {aL,{r) + a^, (r)e-*'"^*) xl, (25) 

where xl is a constant spinor which satisfies j^Xl = —Xl ■ We have tried to be cautious when 
proposing this mixed state, since we are assuming that there is a rest reference frame for the 
particle in such a mixed state. With this reasoning, the Dirac equation i'y^d^ipL,mix = ^n'ipR 
is satisfied, since mi7°XLe~*™^* = miXRC"*™^* {i^Xr = Xr) and there is no right-handed 
massless state, neither inside nor outside the branes. The probability density associated with 
this mixed state is given by 

p{r,t) = \N\^ [aL,{ry + aL,{ry + 2aL,{r)aL ,{r) cos{m^t)] , (26) 

which is an oscillating probability density with period of oscillation T = h/mic^, such that a 
tiny mass implies a long tunneling time. In this scenario the fermion tunnels from one brane 
to the other, being likely found on both branes, but not simultaneously. As L increases, other 
massive localized states can be realized inside the branes. In fact, we have found numerically 
that there is room to one more massive state in the present case. The appearance of a tower 
of locahzed massive states is very dependent on the deepness and width of the double-well 
effective potential UL{r), that is ultimately dependent on the superpotential W{(f)) whose 
classical solution is a kink that can be continuously deformed into two-kink solution . In the 
next section we comment on the construction of such models. 

5 Conclusions 

We have studied the mechanism that leads to the localization of massless fermions on a 
split brane in the cases of warped and flat geometries. The brane is immersed in a five- 
dimensional space-time and is defined by the behavior of a scalar field coupled with gravity 
in the case of warped space-time. The nonpolynomial potential of the self-interacting scalar 
field which generates the split brane was introduced before in reference [18], but any other 
model which has deformable solitons as minimal energy configurations could be used as well, 
for example, a convenient 0^ polynomial potential [25] 

The case of fiat geometry is more manageable than the case of warped geometry, not only 
because we obtain the localized modes in a simple way, but also, and mainly, because it 
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allows us to understand why the coupling of fermions with the scalar field should be chosen 
in such a way that a supersymmetry in a fundamental level is realized. The most convenient 
coupling, which preserves the supersymmetry, also leads to the expected behavior of the 
fermion wavef unctions, that is, the ground-state wavefunction follows the brane splitting. 
As the brane splits into two branes the wavefunction is also split, with peaks at the cores 
of the two branes. Although the simplest Yukawa coupling of fermions to scalar fields also 
provides localization of massless fermions, the corresponding massless fermion wavefunction 
present a peak just between the two branes, such that the observation of massless fermions 
is suppressed inside the branes. Hence, one concludes the supersymmetry inspired coupling 
as the most adequate one to study the entrapment of fermions by a split brane. 

We also note that eventual massive localized states are more difficult to be found in 
the case of warped geometry than in the case of fiat geometry. In the later case we have 
also discussed the tunneling of a massive fermion between the branes, which lasts until 
the branes are infinitely separated from each other. That critical limit is described by 
another model, obtained from the starting one by taking the limit b — ?■ —a in ([6]), that is 
V^cj)) = 2fi^{(f)'^ — a\(f)\y, whose solution is one of the single kinks (j) = |[±H-tanh(/aa(r=F'"o))]! 
with To a reference point where the core of the defect is localized. In this case we found that 
the energy gap between the first localized states is large enough such that they can not be 
considered as quasi-degenerate. It can also be observed that the number of massive localized 
states depends on the deepness and width of the quantum mechanics effective potentials, 
which are defined by the field theory model one chooses to deal with. We have mentioned 
above that others models whose classical solution exhibit a two-kink profile (split brane) can 
afford a tower of localized massive fermions states. A class of such models, called deformed 
models, have been proposed ^7\, as deformation of others known models. We notice that 
those models can also be constructed, together with new ones, from the deformation of 
zero modes excitations of well know models whose classical solutions are single kinks. This 
proposal is being analyzed in more detail and will be reported later elsewhere. 
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R(r) 





Figure 1: Ricci scalar for L = 0.01 (dashed line), L = 1.6 (thin solid line) and L — 4.5 (thick 
solid line) evidences the formation of a double-wall structure as L increases. 
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0.8 r 




Figure 3: ai oi^") (A^-t space-time) in the case F{(p) = —W^^, for L = 1.6 (thin sohd hne) 
and L = 4.5 (thick sohd hne). 
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Figure 4: Effective potentials of equations (123|1 with L = 1.5. ^/^(r) (sohde hne), UR{r) 
(dashed hne). 
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Figure 5: «Li(r) (upper) and «_Ri(r) (lower) in the case -F(0) = —W,f,^, for L = 1.6 (thin 
sohd hne) and L = 4.5 (thick sohd hne), in flat space-time 
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